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?-H ■ The effect of an alternating potential on a one-dimensional half-filled Hubbard model with 

'^ . repulsive interaction has been examined by applying the renormalization group method to 

C^ , the bosonized Hamiltonian. The electronic state, which is determined by the competition 

between alternating potential and umklapp scattering, is calculated where the relevance and 
i-rt \ the irrelevance of the alternating potential leads to the band insulator and the Mott insulator 

g : respectively. The excitation gaps for charge and spin fluctuations are calculated for both states. 
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§1. Introduction 



Q> \ An insulator with mixed ionic-covalent character in the presence of electron correlation is of 

interest since the large repulsive interaction leads to a Mott insulator due to umklapp scattering 
^ ■ at half-filling. 

The effect of electron correlation on such an insulator has been studied in some materials which 

^ \ consist of two-kinds of atoms per cell. In the perovskite oxide compounds, BaTiOa, where the 

Q \ ferroelectricity is usually discussed based on the lattice distortion, the covalency plays an important 

role due to the following fact. The first principle calculation shows that Ba^^Ti^'^^+03' ^ and then 

^ ■ the valence is reduced significantly from the formal valence.EI^ Such a prediction is consistent with 

C^ ' the XPS experiment. a* Further the effect of correlation is also claimed by the fact that dynamical 

effective charges for oxygen and titanium are large due to the covalence.H'o' 

In terms of exact diagonalization, the effect of the electron correlation has been examined for a 
one-dimensional half-filled Hubbard model with alternating potential, which is caused by two kinds 
of atoms and induces the difference of the site energy.H'EP When the repulsive interaction increases 
or the level separation decreases, the band (ionic-covalent) insulator undergoes a transition into 
the Mott insulator followed by a jump of electron numbers of each lattice sites. In addition to 



the electron-lattice interaction, the electronic polar izability is strongly enhanced near the phase 
boundary due to the correlation.LP 

Polarization induced by a sublattice displacement, has been examined for a one-dimensional 
two-band Hubbard model at half-filling. QH-* The number of parameters is reduced by setting equal 
values of the Hubbard U on two kinds of atoms of oxygen and generic cation, while the site energies 
are different. At a critical value of interaction corresponding to a transition from a band insulator 
to a Mott insulator, the static ionic charge is continuous but the polarization and the dynamic 
charge are discontinuous showing the Berry's phased^ associated with the macroscopic polarization 
as the primary order parameter. The discontinuity disappears by removing the finite-size effect, 
i.e., displacing the coarse mesh in order to avoid the A; = singular point and taking the large 
number of A;-points where k is the quasicell momenta. 

The model which is the same as ref. Q but has the infinite length, has been explored by use 
of mean-field approximation in order to clarify the nature of electronic states and the relation of 
the Berry's phase to the phase transition.t3* The ground state has antiferromagnetic ordering for 
large repulsive interaction indicating a magnetic instability while the state is nonmagnetic for small 
interaction. There exists the gap to the charge excitations independently of the value of repulsive 
interaction while the system is replaced by antiferromagnet with gapless spin excitations in strong 
coupling regime and has a gap to all excitations in the opposite limit. Although the magnetic state 
with the gapless spin excitation is known for the conventional Hubbard model with the repulsive 
interaction, it is not clear if such an excitation still exists even for the small limit of the alternating 
potential. Therefore the method beyond the mean-field theory is needed to study the model with 
the infinite length. 

In order to treat strictly one-dimensional fluctuations and electron correlation, we apply, in the 
present paper, a method of the bosonization to the system, which is the same as refs. § and 



lOj . Further, the renormalization group method is utilized where the electronic states of the band 
insulator and the Mott insulator are determined by the relevance and the irrelevance of alternating 
potential. In §2, formulation for deriving the renormalization group equations is given in terms of 
the phase variable, which is based on the bosonization. In §3, the transition from the band insulator 
to the Mott insulator is calculated as the function of on-site electron-electron interaction and the 
energy difference between two-sites. In §4, summary is given and the effect of the alternation of 
the electron hopping energy given by the Su-Schrieffer-Heeger model,EJ^ is briefly discussed. 

§2. Formulation 

The Hamiltonian of one-dimensional half-filled Hubbard model with the alternating potential is 
given by 

^ = -* E (s> '^j+i.- + '^•^•) + ^0 E (-1)' %■- + ^E %-T ^u ' (1) 



where cj^ denotes a creation operator for the electron at the j-th site with spin cr(=t, |). rij^cr = 
c\ „Cj a- Quantities t, Wq and U are energies for the transfer integral, the alternating potential and 
the on-site electron-electron Coulomb repulsive interaction, respectively. 

First, eq. (1) is rewritten by use of the Fourier transform, Ck^a = l/V^V ^- e~^ ^ Cj^o-, where Rj 
is the location of the j'-th lattice site and A^ is the total number of the lattice site. The first and 
second terms of eq. (1) are expressed as, Y.k,a ^k 4,^ Cfc,<7 + Wq Y.k>o,a (4,aCfc-7r/a,<7 + cl-7T/a,a'^k,a) 
where e^ = —2t cos ka and a is the lattice constant. The dispersion of the kinetic energy, e^, around 
the Fermi momentum is linearized as v-pipk — k^), where v-f{= 2ta) and A;f(= it /2a) denote the 
Fermi velocity and the Fermi momentum, respectively. The fermion operator around the Fermi 
momentum is rewritten as Ck,p,cr{= Ck+pkp,^), where p{= + and — ) denotes the branch for the right 
moving and left moving electrons. The third term of eq. (1) is rewritten as (L = Na) 

-xy _ 51 Y^ t t 

ki,k2,q,P 

, 52 v^ t t 

"I" ^ Z^ '^ki,p,^'^k2-p,i^f'2+q,-p,l'^ki-q,p,'i 
ki,k2,g,P 

, S3 sr^ \ t 

'^ j^ Z^ '^ki,p,l;'^k2,p,l'^k2+q,-p,i^ki-q,-p,'i 
ki,k2,g,P 

+ ^ Z^ '^ki,p,\'^k2,p,i'^k2+q,p,i'^ki~q,p,t , (2) 

ki,k2,q,P 

where gi = §2 = 93 = 94 = Ua and these quantities denote coupling constants for the backward 
scattering {91), forward scattering between the opposite branch (5^2), umklapp scattering (33) and 
forward scattering within the same branch (54), respectively. 

Next the bosonization methodtJ^ is applied to eq. (1). The phase variables, dp± and 9a±, ex- 
pressing fluctuations of the charge density and spin density are given by 

^p±(^) = 4f E ^^"^'"^'"''^ E (4+,,+,.cfc,+,<. ± 4+,,-,.^^,-,^ > (3) 

^ q^O k,a 

^<^±(^) = 4 E i^"^'"'"''^ E ^ (4+,,+,.c^,+,<. ± i-,q,-,.'^k,-,.) , (4) 

V^^^O^ k,a 

where a = +(— ) corresponds to t(|) and [6u+{x),9,^i^{x')] = mS^y sgn{x — x'). By using these 
phase variables, we rewrite the field operator, ■0p^o-(a;)(= L^^ J2k ^^^''^^'^^^^ '^k,p,a) , as 

'^p<" = ~7^= ^^P i^pf^FX + iQp.cr) exp (i-TrHp,^) , (5) 

\/2'Ka 

Qp,a = ^[p0p++9p^+a{pe^++e^-)] , (6) 

where a is of the order of the lattice constant. The phase factor in eq. (5), is chosen as Hp^+ = 
p(iV+,+ + N-,+)/2 and Hp _ = {N+,+ + N-,+) + p{N+- + N--)/2 with the number operator, 
Np^fj, in order to satisfy the anticommutation relation. Such a choice of Hp^o- conserves a sign of the 



nonlinear terms, which is obtained from the interaction terms of eq. (2). Based on these bosonic 
fields, eq. (1) is rewritten agl3 
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of eq. (7). From the assumption that response functions are invariant for a — > q' 
renormalization group equations for Kp, K(^,g3,gi and Wq are derived as (Appendix A) 
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dl 



where Vp = [{Ittvf + 94? - glV^^f^T^, Kp = [(27r?;F + 54 - 52)/(2^fF +94+ 52)]^/^ v 
94? - 5?]^/V2vr and K, = [(27r7;F - 54 + 5i)/(2vrfF - 54 - 5i)]^/^- 

Now we derive renormalization group equations by examining response functions, Ra{xi —X2,ti — 
T2) = {TtOa{xi,ti) 0\{x2,T2)) , where Tj is the imaginary time and Oa denotes several kinds of 
order parameters. Response functions are evaluated perturbatively by expanding nonlinear terms 

(8) 
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where initial conditions are chosen as Ki,{Q) = K^ {u = p and a), G3(0) = g^/2TTv-p, 6*1(0) = 
gi/2iTVY and Gvk(O) = Wo/{vYa~^). For the simplicity, we replaced Vp and v^j by wf, in eqs. (8)- 
(12) since the deviation of the numerical factor is small. 

The calculation of K„{1) in eq. (9) is performed by the following procedure. When we substitute 
K^{1) = 1 + Gi{l) for r.h.s. of eqs. (9) and (11) and retain the coupling constants up to the second 
order, eq. (9) becomes equal to eq. (11). In this case, the fixed point of K^{1) with Wq = is given 
by unity showing the validity of the SU(2) symmetry, i.e, Kfj[oo) = 1. Therefore eqs. (9) and (11) 
are treated by such an expanded form for K„[l) > 1 while we use the original form of eq. (9) for 
K^il) < 1. 

We note the electronic states obtained from eq. (12). In the absence of the [/-term, the state of 
eq. (1) has a gap, 2Wq, at k = +kF of the band while that of eq. (1) with Wq = exhibits a charge 
gap, Ap, due to half-filling corresponding to a Mott gap. Both of these cases are insulating states. 
The former is the band insulator with the relevant Gwil), which increases infinity for large /. The 



latter is the Mott insulator with the relevant G3{1). In the presence of both U-tema and Wo'term, 
there are two cases of the relevant Gwi^) and the irrelevant Gwi^)- We use the terminology of the 
band insulator for the relevant Gwil) and the Mott insulator for the irrelevant Gwil) although 
the latter case denotes actually a crossover to the Mott state, i.e, the spin gap decreases rapidly to 
zero as shown in the next section. 

Here we examine the order parameters, which are characteristic of these insulating states. Since 
the phase of 6p^ is locked at vr or (±7r/2) for the relevant (irrelevant) Gvi/(0) '^^ can consider 
the order parameter of the band insulator (the Mott insulator) as (cos0p+) ((sin^p^)). The finite 
value of (cos 0p+) denotes the ionic state since the Wo-term of eqs. (1) and (7) is proportional to 
the difference between electron numbers of two kinds of lattice sites. As for the locking of ^ct+7 
the effect of the Wo-term is much stronger than that of the ^fi-term in eq. (7), because the g(i-term 
becomes irrelevant for Wq = 0. Thus we always have the finite value of (cos0(j+)- We note that 
the quantity (cos 6'p+) {cos9a+) is negative due to the positive value of Wq in eq. (7) and that 
there is a degeneracy of a state with the positive (cos0p+) and the negative (cos0(j+) and a state 
with the negative (cos0p+) and the positive (cos 00-+)- For studying these properties, we calculate 
response functions, Rcosd^+ifi — ^2) = (TT-cos0;y+(ri) cos0i/+(r2)) {v = p and a). By use of the 
renormalization group technique, these response functions are derived as (Appendix A) 

Anir/a) 

/o 

An{r/a) 

/o 

where r = [x^ + (wft)^]^'^. For the band insulator the limiting values of both RcosO +(^) and 
Rcose^+i'>') are finite, while only that of Rcose^+i'>') becomes finite for the Mott insulator. 

§3. Band insulator vs. Mott insulator 

For the calculation of characteristic energy in terms of renormalization group equations, it is 
crucial to determine the magnitude of the cutoff parameter, a, which is chosen conventionally as 
a = a/vr.llj) We estimate a from the response function for charge density wave with a momentum 
close to 2A;f as follows. Such a response function with Wq = U = can be calculated explicitly by 
both the conventional method and the bosonization method (Appendix B). From the comparison of 
these two quantities for the half-filled band, we obtain a ~ 7r/(1.273a). Thus characteristic energy 
corresponding to / is given by vpa^^ exp[— /]. 

In Fig. 1, quantities i^p(/)(solid curve), K(j{l) {dashed curve), G3(/) (dot-dashed curve) and 
Gpi/(/) (dotted curve) are shown as a function of / for U/t = 2 and Wo/t = 1. With increasing 
/, Gwil) increases monotonically and then becomes relevant leading to the band insulator. The 
change of sign of G3(/) occurs at / ~ 2.1 due to the relevant Giy(/). The positive G^i^l) originates in 
the umklapp scattering, 5^3, and the negative G^il) comes from the higher harmonics induced by the 



Rcosep+ir) =exp 
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' dl (Kp(l) + 2G^{1) 
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Fig. 1. The /-dependence of coupling constants, Kp{l), Kcr{l), G-i{l) and Gwil), with fixed U/t = 2 and Wo/t = 1. 
The 2-dot-dashed curve denotes Kp{l) for Wo/t = 0. 
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Fig. 2. The /-dependence of coupling constants, Kp{l), K„{1), G3{1) and Gw{l), with fixed U/t — 4 and Wo/t — 1 
where the notations are the same as Fig. 1. 



VFo-term. With increasing I, Kp decreases monotonically since both G'i{l) and Gw{l) reduce Kp{l). 
The reduction of Kp{l) by Gw{l) is understood from the fact that Kp{l) is suppressed noticeably 
compared with the case of Wq = (2-dot-dashed curve). The /-dependence of Gi{l) is similar to 
that of Gz{l)-, i.e., the positive Gi{l) changes the sign with increasing /. With increasing /, Kp[l) 
and K(j[l) merge each other and their magnitudes are reduced to zero. This behavior indicates the 
formation of both charge gap, Eg, and spin gap, Ao-, which originate in the VFo-term. 

In Fig. 2, the numerical result for U /t = 4 and W^/t = 1 is shown where the notations are the 
same as Fig. 1. The behavior of K„{1), G3(/) and G\y{l) is quite different compared with that of 
Fig. 1 although that of Kp{l) is qualitatively the same. The difference between Kp{l) (solid curve) 
and that of Wq = (2-dot-dashed curve) is very small due to the irrelevant Gw{l)- The effect of 
the Wo-term. reducing the charge gap is seen from the fact that Kp{l) in the presence of Wq (solid 
curve) is slightly larger than that of Wq = (2-dot-dashed curve). With increasing /, Gw{l) takes 
a maximum and decreases to zero. Thus the effect of Wq becomes negligibly small for the large I 
(i.e., small energy). Such an irrelevant Gw{^) leads to a state similar to the Mott insulator induced 
by the umklapp scattering. The quantity G3(/) is always positive and increases monotonically from 
the initial value 53. The /-dependence of Ka-{1) is similar to that with Wq = 0, i.e., Ka-{1) ^ 1 in 
the range of Fig. 2. However a small spin gap does exist since Ku{l) reduces to zero for the large /. 
Based on the calculation given by Figs. 1 and 2, it turns out that there is a critical value ofU = Uc, 
where the band insulator is obtained for U < Uc and the Mott insulator is obtained for U > Uc- 
In Fig. 3, the fZ-dependence of excitation gap and | {cos6u+) \ is shown. The charge gap (Eg) and 
spin gap (Act) are evaluated by the formula. Eg = vpa'^ exp[— /p] and Act = vpa^^ exp[— /ct] where 
ly (u = p and a) is chosen as Ky{ly) = 0.3. Such a choice is reasonable since, in the absence of 
Wq = 0, the charge gap of the present calculation (dotted curve) coincides well with the exact one 
(dot-dashed curve) .tj) The [/-dependence of Eg and Act is shown by solid curve and dashed curve, 
respectively. The arrow denotes the location for a critical value of Uc/t ~ 3.0, which separates 
the band insulator from the Mott insulator. There is a minimum of Eg (solid curve) where the 
corresponding U is slightly smaller than Uc- Such a minimum arises from a competition between 
the PVo-term and the (73-term. For large U , the magnitude of Eg moves close to Ap(dotted curve), 
which denotes the charge gap in the absence of Wq. With increasing [/, the spin gap Act decreases 
but takes a finite value even for U > Uc- For small [/, one finds Eg ~ Act indicating that the 
difference between the magnitude of the charge gap and that of the spin gap is negligibly small. 
Thus it turns out that the behavior of the charge gap is quite different from that of spin gap. 

In the inset of Fig. 3, the [/-dependence of | (cos 0p+) | and | (cos 9^+) \ is shown by solid curve (1) 
and dashed curve (2), respectively where the arrow denotes U = Uc- The quantity (cos0;^+) {v = p 
and o"), is defined by (TT-cos6'jy+(x) cos0,y+(O)) -^ (cos0i/+) for |x| -^ 00. In the context of the 
present calculation, such a value is calculated from the minimum value of response function due 




U It 

Fig. 3. The [/-dependence of Eg/t and Ao-/i is shown with fixed Wo/t — 1. The arrow denotes a location for a 
critical value, Uc/t{^ 3.0), where the band insulator is obtained for U < Uc and the Mott insulator is obtained 
for U > Uc- The charge gap, Ap/t, for Wo/t = is shown by the dotted curve (the present calculation) and the 
dot-dashed curve (the exact solution) 113 In the inset, the f/-dependence of | (cos 6'p+) | and | {cos 6a-+) j is shown 
by curve (1) and curve (2) respectively where the arrow corresponds to Uc- 



to the second order renormalization group equation. With increasing U/t, \ (cos 6p-\-) \ decreases 
and is reduced to zero at U = Uc, corresponding to the transition from the band insulator to the 
Mott insulator, where the locking of phase for the former (latter) state is given by 9p^ = vr or 
{9p-\- = ±it/2)- The finite value of (cos0p+) corresponds to the relevance of Gh/(/), i.e, the band 
insulator. In contrast to (cos 0p+) , we find that (sin 9p^) becomes finite for the irrelevance of Gw{l), 
i.e., U > Uc due to the change of the locked 6p+ from to ±-k/2. The quantity | (cos0ct+) | decreases 
but remains finite even for U > Uc showing the similarity to A„. We note that | (cos ^p+) | oc ^|Wo| 
for [7 = and small |VFo|/t ^ 2. Such a result is compatible with the fact that the conventional band 
gap with [7 = 0, is proportional to Wq. Then the present calculation is valid for Wo/t ^2. The 
[/-dependence of Eg and (cos0p+) (Aq- and (cos0ct+)) shows that a transition (a crossover) from 
the band insulator to the Mott insulator appears due to the competition between the alternating 
potential and the electron-electron interaction. 

From the calculation of Uc similar to Fig. 3 with several choices of Wq, a phase diagram of the 
band insulator and Mott insulator is shown on the plane of U/t and Wo/t in Fig. 4. Since the the 
boundary between these two states is convex downward, the effect of Wo is much larger than U 
for small U. The boundary, which is well reproduced by a formula, exp[—2Trt/U], for small U/t, 
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Fig. 4. Phase diagram for the band insulator and the Mott insulator on the plane of C//f and Wo/t. The inset shows 
the same phase diagram on the plane of Ap/t and Wo/t. 



indicates a fact that the characteristic energy for the competition is not U but the charge gap Ap. 
In the inset, the same phase diagram is shown on the plane of Ap/t and Wo/t. From the result 
that Wo/Ap = 1.6 ~ 1.7 for 0.3 < Ap/t < 2.0, it is found that the boundary between the band 
insulator and Mott insulator is determined by the competition between Ap and Wq. 

§4. Summary and Discussion 

In the present paper, we have examined the effect of the alternating potential, Wq, on a one- 
dimensional half-filled Hubbard model with the repulsive interaction, U. The property of the model 
with the infinite length is calculated by use of the methods of bosonization and renormalization 
group. With increasing U, a second order phase transition followed by the continuous change of the 
order parameters occurs from the band insulator into the Mott insulator at Uc where the former 
(latter) corresponds to the relevant (irrelevant) Gwi^)- Such a transition is demonstrated from the 
fact that (cos 0p+) becomes zero for U ^ Uc — 0. The band insulator denotes the ionic state since 
(cos^p+) 7^ in addition to (cos0o-+) 7^ 0. Phases of the band insulator are locked at 9p^ = vr 
(or 0) and ^0-+ = (or tt) respectively. For the Mott insulator, phases are locked at 9p^ = tt/2 or 
—it/2 and 0o-+ = or vr due to the presence of the spin gap even if Gwi^) is irrelevant. The charge 
gap. Eg, takes a minimum and the spin gap A^- exhibits a crossover around the transition point of 
U = Uc- The phase diagram of band insulator and Mott insulator has been obtained on the plane 
of U and Wq. 



Here we compare the boundary calculated in the present paper with those obtained by the 
numerical diagonalization. In ref. ^, in which the same model is treated, the boundary between 
band insulator and Mott insulator is calculated as Uc/t ~ 2.27 for Wo/t ~ 0.57 while our result 
shows Uc/t ~ 2.4 for W^/t = 0.57. In refs. |5| and ^ which introduced two kinds on-site potentials 
Ub < Ua in addition to the present model, they obtained Usc/t — 2.5, 3.4 and 5.4 for Wo/t ~ 
0.67, 1.0 and 2.0 while the present results show Uc/t ~ 2.5, 3.0 and 4.2 for Wo/t = 0.67, 1.0 and 2.0 
respectively. A quantitative coincidence between their numerical calculations and ours are found 
for Wo/t ^ 1. Further we note a fact deduced from the present result that (cos 6*^+) = at C/ = Uc- 
Close to U = Uc, the free energy as a function of (cos 0p+) takes a double minimum for the band 
insulator due to the degeneracy with respect to 9p^ = and vr while it takes a single minimum for 
the Mott insulator. Thus the charge fluctuation is strongly enhanced around U = Uc as seen in the 
second order phase transition. Such a fact could be related to the enhancement of the dynamical 
charges' and the increase of localization lengthH' at the critical point. 

Finally we discuss about the effect of dimerization where the corresponding Hamiltonian is given 
by 

-tdE [(-!)' sWi,.h.c.] . (15) 

In terms of phase variable, eqs. (5) and (6), one can replace eq. (15) as 

4td 



dx sin 0p+ cos 6*0-+ . (16) 

■Ka J 

Equation (16) is compatible with the (73-term in eq. (7) since the energy gain of the former is 
obtained for 0p+ = 7r/2 and ^0-+ = vr (or 6pj^ = — 7r/2 and 6^+ = 0) while that of the latter is 
obtained for 9pj^ = ib7r/2. By noting that td-term and the PVo-term have the same periodicity with 
27r as a function of 9p+, these two terms coexist for the relevant Gw{l)- In this case, the phase is 
determined so as to minimize two energies of t^-iexva. and Wo-tsrm. We can incorporate eq. (16) 
into renormalization equations and find that td-term is always relevant and that there exists a 
boundary between the relevant Gw{l) and the irrelevant Gw{l)- Based on these consideration, we 
can expect following. Although the qualitative feature is the same as that with t^ = 0, there are 
noticeable facts that the spin gap is strongly enhanced due to the relevant t^ and that the critical 
value Uc decreases due to the enhancement of umklapp scattering by td- 
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Appendix A: Derivation of the renormalization group equations 

A.l Coupling Constants 

The renormalization group equations are derived in a way similar to Giamarchi and Schulz.EJ) 
By treating the nonlinear terms in eq. (7) as the perturbation, the response function, Rp{x,T) = 
(TT-exp[i0p+(x,r)] exp[— i0p+(O, 0)]), is calculated up to the third order as 
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'0 \ /OJ \ /o 

IGsG^^v- /• dVs dV dVs 



e,e 



gi[ep+{ri)-ep+(r2)]\ /gie[2ep+(r3)-ep+(r4)-ep+(r5)]\ I /gie'[e,+ (r4)-0<.+ (r5) 



' , J o? o? a2 

i[ep+{ri)-ep+(r2)]\ / ie[ep+{r4)-ep+{r5)]\ \ / i€'[2e,+ (r3)-e.+ (r4)-e<.+ (r5 



'W. V /" 1_[3 !llilZl r / i[0p+{ri)-ep+{r2)+eep+(r4)-eep+(r5)] 

2 (2^)3 f^^y a2 ^^2 ^2 W /o 



+ ••• , (A-l) 



\ /OJ \ /o 



where r = (xjVft) and (• • •)o denotes an average over the harmonic parts including Kp or K^, in 
eq. (7). Scaling equations for coupling constants up to the second-order can be calculated from 
response function, Rp(r), which is expanded up to the third-order for the nonlinear terms. A 
straightforward calculation of eq. (A-l) yields. 



Rp{ri - r2) = e 



-KpU(n-r2) 



1 Gi ^ fd^rgdV _K 



o? o? 



e 



' pU(TY-r-i)-^KpV{:rz-T^) 



^ 2 (2vrV ,2, 

^ rg2eKp[C/(ri-r3)-C/(ri-r4)-C/(r2-r3)+C/(r2-r4)] _ A 

e 

^ rg2eXp[t/(ri-r3)-C/(r2-r3)]+ei^p[-C/(ri-r4)-i7(ri-r5)+(7(r2-r4)+(7(r2-r5)] _ -A 



^W V^ /"^L^^Lilp--f^pC/{ri-r2)-(i^p+Ka)C/{r3-r4) 
{27t)^ £^J a2 a2 



GsC^y -^ /■ d rsd r4d r5 ^_;^^^(^^_^^)_^(j^^_j^_^)^(^^_^^)^_2j:fjj7(r3-r.4)+L/(r3-r5) 
(2^)3 ^J a2 „2 ^2 
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(27r)3 ^y a2 ^2 ^2 « 

X fQ^Kp[U{ri~r4)-Uiri-r5)-U{r2-r4)+U{r2-r5)] _ A ^ ^ (A-2) 

where U{r) = ln(r/a) (r = [x^ + (vft)'^]^''^). We have discarded the difference of velocities, i.e., 
Vp = Va = vp. As for third-order with respect to couphng constant in eq. (A- 2), three cases of 
rs = r4 + r, rs = r3 + r and r4 = ra + r with small r are chosen where the expansion of these terms 
around r = is rewritten as 

_ GsGy^^ ^ r d r3 d r4 K^u(^ri-r2)-4KpU{r;i-r4) 

(27r)3 ^J a2 c^2 « 

^ rg2ei^p[(7{ri-r3)-;7(ri-r4)-;7{r2-r3)+C/{r2-r4)] _ A /" ^ ?" ^(K,-K^)U(r) 
_ n ^sGw V- /■ ^^^^ -ii'p!7(ri-r2)-{i<'p+K^)l/(r3-r4) 

(27r)3 ^7 a2 „2 « 

X ff,'^Kp[U{ri-r:i)-U{ri-ri)-Uir2-r:i)+U{r2-ri)] _ ^\ f ^ r ^_2K„U(r) 



X L^Kp[U(ri-r3)-U(ri-r4)-~U(r2~r3)+U(r2-r4)] _ A f ^ ^-2K^U{r) + . . . _ (A-3) 



(27r)3 ^J a2 ^2 

d^r 

By comparing eqs. (A-2) with (A-3), it is found that eq. (A-3) has an effect of the renormalization 
for Gl-term and G^-term in eq. (A-2). Actually eq. (A-2) can be rewritten by the first three terms 
with effective coupling constants, which are given by 



°5'-°3-Oi-/$(£)'^''-'-'\ ^ (A.4) 

Gw = Gw — Gs Gw I — 2" ( ~ ) - Gi Gw / — 2" ( ~ ) • (^'5) 

On the other hand, the effective quantity of Kp can be calculated from the response function up 
to the second order for the nonlinear terms. By expanding the exponential in the second and the 
third terms of eq. (A-2) and using new variables r = rs — r4 and R = {r^ + r4)/2, eq. (A-2) is 
rewritten as, 



Rpin -r2) = e 



-KpU(ri-r2) 



2 



(2vr) 
where 



l + 2^K'^J+{2Kp)mn-r2) 

+ 7&^P J+(.Kp/2 + KJ2) /+(ri - r^] 



(A-6) 



I+{ri -r2) = I d?R U{ri - R) V| U{r2 - R) = 27t U{n - r2) , (A-7) 

^.(i^)= /^4e-^"^^^) = 2. r^f^)"" . (A-8) 
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The scale invariance of eq. (A-6) is shown by reexponentiating eq. (A-6) where the effective quantity 



of Kp is obtained as 



Kf = Kp-2GlK'J^(^) "''''' -Gl.K'J^(^) "^'"''''^ , (A- 



as 



For the transformation given by a — > a' = ae*^',Ej) these quantities have to be scaled 



Kf{Kl,G';a') = Kf{K,,G^-a) , (A-10) 






Gf{K,G'^-a') = Gf{K,,G,-a){a'/aY-'^^ , (A-11) 

Gl^{K,G'^-a') = G'^{K,,G^;a) {a' /af-""^/^-^^'^ , (A-12) 

where K^v = p, a) and G' {p = 3, 1, W) denote renormalized quantities. The quantity 2Kp {Kp/2+ 
Kfj/2) is a scaling dimension for G'^ (for C^), which is obtained from the second (third) term 
of eq. (A-9). By applying the infinitesimal transform to eqs. (A-9), (A-4) and (A-5), renormalized 
quantities are expressed as 

K'p = Kp-2GlKldl-GlrKldl , (A- 13) 

G^ = G3 + (2 - 2Kp) G3 d/ - Gl, dl , (A- 14) 

G'w = Gw + {2- Kp/2 - K^/2) Gw dl - G3 Gw dl - Gi Gw dl , (A- 15) 

which lead to eqs. (8), (10) and (12), respectively. The renormalization group equations for K„{1) 
and Gi{l) can be obtained in a similar way to eqs. (8) and (10) by calculating the response function 
given by (T,- exp[i0o-+(2;)''")] exp[— i0o-+(0, 0)]). We note that, in case of Wo/t = 0, these equations 
reduces to to well known equations for half-filled Hubbard model.113) 

A. 2 Response Function 

Based on the solution of of eqs. (8)-(12), we examine the response function, Rcose +(^^1 — ^2) = 
{Tr cos0p+{ri) cos6p+{r2)), which is calculated by writing Rcosep+{r) = ex.p[-KpU{r)] ■ Fcose^+ir)- 
The perturbative calculation for Rcose + leads to 

G3 /-dVg 



^cose,+ [ri r2j-i ^^ J ^^ e 

Kp{n-r2) (A-16) 



Ht" / T" \ ^^p f c\t' / T" \ P ^cr 



2Gl -(-) +GI, 

The scaling invariance under the transformation of a ^ a' = oe is assumed for the quantity 
^cos ep+ (r) , which is expressed as 

Fcos9,+ {r,K,,G^;a) = hos9,+ {dl,K,,G^) ■ Fcos9,+ {^,K,G'^;a') . (A-17) 

The quantity /cos0p+ is determined perturbatively. By applying the infinitesimal transform to 
eq. (AT6), we obtain 

Fcos e,+ = 1 - 2 G3 d/ + (2 Gl + Gl,) K^ U{n - ^2) dl 
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G', f dVs 



2tt J a' a'2 



,2i^p[(7'{ri-r2)-C/'{ri-r3)-C/'(r2-r3)] 



OO A^ / /r \ 3—4:K, 

^ ' .' a' \a' 






KlU'{ri-r2) ,(A-18) 



where U'{r) = ln(r/a'). From eq. (A-18), the multiplicative factor /cose + is obtained as 



/cose,+ (d/, K., G^) = exp -2 G3 d/ + 2 G^ K' U{r) dl + G^ /C^ C/(r) d/ 



(A.19) 



The transformation of eq. (A- 17) is performed recursively until the new cutoff a' reaches r. Thus 
the reconstructed Fcosep+ is expressed as 

f\n(r/a) 



Fcos 0p+ {r,K„,Gf,;a) = exp 



dnnI,,,g^4dl,K,{l),G^{l)) 



(A-20) 



The terms including the second order of the coupling constants are rewritten as follows. From 
eq. (8), one obtains, 

r>ln(r/o) 



/ d/ 

Jo 

= Kp{0)\n{r/a) 



2Gi{l)K'p{l) + Giy{l)K'p{l)\\n 

/•ln(r/a) 

/ d/Kp(/), 

/o 

where a{l) = aeK Thus the response function -Rcos6»p+ is expressed as 

fln(r/o) 



a{l) 



(A-21) 



Rcosep+{r) = exp 



(A-22) 



j dl(Kp{l) + 2G^{l)] 

which leads to eq. (13). In a similar way, another response function leads to eqs. (14). 

Appendix B: Evaluation of cutoff parameter 

We calculate response function of charge density wave with Q = 2kp + q and kp = it /(2a), which 
is expressed as 

n(Q) = \J2J^^ (r4^(T)c,,.(T) 4,^(0) co,.(0))e- 



'iQRj+iuJmT 



J.o- 



iuJm^O 



(B-l) 



For the case with only the first term of eq. (1), eq. (BT) is calculated as 



n(Q) = J-r dkl^^^^±^l^I^ -^ : 

27r J-n/a Sk - Ek+Q ga<i 2'KVY qa 



' i„^ 



(B.2) 



where C ~ 1.2737r. The evaluation by the bosonization method leads t^^^' 

1 /■ . , 1 



nB(Q) 



(27ra)2 
Ko{qa) 



dxdr 



1 + (vpT/a) + (x/a) 



^iqX-lUmT 



iulrn-^O 



' Ini- 



27rvF 9^0 27rt'F qa 
Prom eqs. (B-2) and (B-3), it is found that a ~ a/(1.2737r). 



(B-3) 
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